Electromagnetically induced left-handedness in a dense gas of three level atoms 
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We discuss how a three level system can be used to change the frequency dependent magnetic 
permeability of an atomic gas to be significantly different from one. We derive the conditions for 
such a scheme to be successful and briefly discuss the resulting macroscopic electrodynamics. We 
find that it may be possible to obtain left handed electrodynamics for an atomic gas using three 
atomic levels. 
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I. INTRODUCTION 

Changing the propagation properties of light by de- 
signing a novel material is of interest both from a ba- 
sic science point of view, and for technological applica- 
tions. Recent advances such as slowing downQ or stop- 
ping light 121, or left handed metamaterials Q promise 
advances in fields ranging from optics [J| to quantum 
computation|5j. It is desirable to find new materials in 
which electromagnetic waves exhibit novel behavior, and 
there is a flurry of activity both theoretically and exper- 
imentally in this direction. 

The macroscopic electromagnetic (EM) properties of 
a medium are characterized by the frequency dependent 
dielectric constant ^6] 



D = e{uj)E 
and the magnetic permeability 

B = fi{Lu)H. 



(1) 



(2) 



Speed of an EM wave of frequency lu in this medium 
is given by 



(3) 



VmMeH 



where c is the speed of light in vacuum. The index of 
refraction is then n = ^J ^{ijj)e(bj) . 

The dielectric constant of the medium shows large vari- 
ations near a resonance, i.e. when the frequency of the 
external field is near an internal state transition. This 
makes it easy to change the refractive index of a medium 
by properly tuning the frequency of the EM wave to just 
above or below a transition. Recently this fact was em- 
ployed combined with quantum coherent effects to get 
very high refractive indices in atomic gases 0, Q • 
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Although e{ijj) can change appreciably for a gas of 
atoms at optical frequencies, the magnetic permeability 
/.t(w) is always very close to its free space value. One 
can give arguments in classical electrodynamics to ex- 
plain this 8], or understand it in terms of atomic tran- 
sitions as follows. Magnetic field component of an EM 
wave couples to the atom much weaker than the electric 
field component. The magnetic coupling to an atom is 
proportional to the Bohr magneton /zs = 2m c ~ Q^c^Oi 
while the electric coupling is eao- The fine structure con- 
stant a ~ 1/137 also shows itself in the induced magnetic 
dipole moment. Overall the effect of an EM wave on mag- 
netic permeability is a^ weaker than its effect on the elec- 
tric susceptibility. Another important fact is that mag- 
netic dipole transitions are allowed only between states 
which have the same radial wavefunction, and generally 
two such states are not separated by optical frequencies 
in energy. 

Now that it is hard to get /^(w) to be different than 
one, we need to question why it is important to have 
another value for it. After all it seems from Eq.lPJ) 
that all the optical properties of the medium depend 
on the product ^[ijj)e{uj). The answer to this question 
is that the refractive index alone does not completely 
represent the medium [3. One can imagine two media, 
one with £1(0;) > 0, /xi(w) > and £2(0;) = —€i{uj), 
H2iijj) = —fJ-ii'^)- They would have the same refrac- 
tion index, however, quite different optical properties. 
Materials with both e < and ^ < 0, are called left 
handed materials, named for parity of the coordinate 
frame formed by {E,H,k}. Optical properties of left 
handed and right handed materials differ mainly because, 
the Poynting vector points opposite to k in left handed 
materials. Most remarkable change happens at the inter- 
face between a left handed material and a right handed 
material, where the usual Snell's refraction law gets a 
sign change. In addition to the inverse Snell's law, the 
reverse Cerenkov radiation and the reverse Doppler shift 
would also be possible in such materials [9|. 

Left handed artificial materials in the microwave re- 
gion have recently been built [i3| by assembling a com- 
posite lattice of metallic split ring resonators and metallic 



wires |ll|, with periodicity much smaller than the wave- 
length of the electromagnetic field, or using anomalous 
propagation properties of light in a photonic crystal, 
with periodicity is in the order of the wavelength of the 
electromagnetic radiation _12|. All such systems (called 
metamaterials), require delicate manufacturing of spa- 
tially periodic structures. In the microwave region, im- 
provements of focusing, filtering and steering properties 
of microwaves would be useful for many practical appli- 
cations. Similar improvements would also be valuable for 
applications operating at optical frequencies. In this pa- 
per we examine the case of an atomic gas without any 
spatial periodicity that could exhibit behavior similar to 
metamaterials at optical frequencies. 

We have remarked that the magnetic dipole response 
to an oscillating magnetic field is smaller by a factor of a^ 
compared to the electric dipole response to an oscillating 
magnetic field. In an EM wave E and B fields are always 
perpendicular to each other and are always in phase. If 
one can get the atom to respond to an electric field E 
with a magnetic moment fl perpendicular to it, and in 
phase with it, one can effectively think that the magnetic 
dipole moment is induced by the magnetic field of the 
EM wave. Thus, it is possible get a magnetic response 
which is only a times smaller than the electric response. 
Such response enables one to achieve a regime where the 
propagation properties of light are significantly different. 

The aim of this paper is to explore the feasibility of this 
idea to modify the magnetic permeability of an atomic 
gas electromagnetically. To this end, we introduce a 
model system in the next section and find the neces- 
sary conditions for the applicability of our scheme. In 
section III we present the results of our calculations for 
two different parameter regimes, a dilute gas and a dense 
gas. We then go on to discuss the consequences of our 
results for experiments. Finally, we give a summary of 
our results and conclusions in section IV. 



II. MODEL SYSTEM 

In this section, we construct a model system for which 
the magnetic permeability can be optically modified. We 
also describe the scheme for modification in detail, and 
discuss its limitations. 

One can readily conclude by parity arguments that it 
is not possible to get a magnetic response to an electric 
field if only two states are involved. An electric field 
causes transitions to states which are of opposite parity 
to ground state, and such states do not have a magnetic 
dipole matrix element with the ground state. To over- 
come this difficulty, we use a three level scheme, similar to 
the one used in electromagnetically induced transparency 
(EIT) [3, [13 where an optically thick substance is made 
transparent and exhibits large dispersive response to the 
external field close to atomic resonance. 

The particular EIT scheme here serves several useful 
features required for left-handedness, such as being dis- 



persive and exhibiting resonance phenomena. EIT ma- 
terials do not suffer from linear absorption at resonance. 
They exhibit small transmission losses even at high densi- 
ties. As a consequence of the resonance, the EIT medium 
stores large amount of energy over the cycles of interac- 
tion, leading to strong material response. In order to 
have a negative electric and magnetic material response, 
we need both the macroscopic polarization and magneti- 
zation of the material become simultaneously so strong 
that they would be immune even the sign changes of the 
applied fields. For the weak probe beam, EIT cannot 
achieve this feat single-handedly. By considering a dense 
medium, with many particles within a cubic resonance 
wavelength, we let the local fields in the substance help 
to enhance the material responses. Indeed, we see that, 
circularly polarized probe electric field, under EIT condi- 
tions, together with the help of the Lorentz-Lorenz local 
field contribution, could maintain strong local currents 
that could give rise to large enough magnetization, in- 
sensitive to sign changes of the probe magnetic field. At 
the same time, the electric response also becomes neg- 
ative. The remainder of the section presents the math- 
ematics behind these ideas as well as the conditions of 
their applicability. 

We require the three states to have the following non- 
zero matrix elements : 



(l|erl3) ^ 
(2|erl3) ^ 

{i\m + 0. 



(4) 



Here /i is the magnetic dipole moment operator given by 



M 



Mb 



qlL + gsS + gil 



(5) 



where the first two terms are the magnetic moments due 
to the electronic orbital angular momentum L and spin 
angular momentum S, while the last term is the contri- 
bution of nuclear spin angular momentum /. The co- 
efficients are g^ = 1 and gs ~ 2 (within a small 0.1% 
correction found by quantum electrodynamical calcula- 
tions). Nucleon magneton is about 1800 times smaller 
than the Bohr magneton. Typical nuclear magnetic mo- 
ments are about 1000 times smaller than their electronic 
counterparts and hence usually negligible. If the Hamil- 
tonian is parity invariant, we can choose all the states to 
be eigenstates of the parity operator V. To satisfy the 
requirement Q) one should have 



(1|7'|1) = (2|7'|2) = -(3|7'|3). 



(6) 



We assume that the states |2) and |3) are coupled with 
an intense coherent beam while a weak probe beam will 
excite transitions between |1) and |3). We will investigate 
the dielectric permittivity and magnetic permeability of 
a medium consisting of such atoms as a response to the 
probe beam. 

Such a system of a three level atom interacting with 
those two optical fields in A scheme as depicted in FigQ] 




FIG. 1: Three-level atom interacting with the probe and the 
coupling fields in A scheme as described in the text. 



is described by a Hamiltonian in the form 

H = Ho+Hi, (7) 

where 



and 



Ho — 2_^ h^iRi: 



Hi--- ^ (r!,e-''''*i?3. + c 



(8) 



(9) 



Here, fiLHi arc the energy levels of a free atom, and 



R^ 



\i){i\ are atomic projection operators. The in- 



teraction Hamiltonian is written under the electric dipole 
approximation. The Rabi frequencies associated with the 
optical transitions are defined by 



n. 



dzi ■ Ei 



(10) 



where the Ei stands for the complex amplitude of the 
positive frequency component electric field of the probe 
laser. The electric dipole operator is expressed as 



d3i = e(3 I r I i). 



(11) 



Within the semiclassical theory of optical interactions, 
density matrix of the system evolves according to the 
Liouville equation 



t->-^-h^-y 



(12) 



Assuming a diagonal relaxation matrix {i \ T \ j) — ^idtj, 
for our model Hamiltonian, density matrix equations (op- 
tical Bloch equations) become |7l[l^ 

P33 = -73P33 - ^ ^ (f^re-'Vsz - c.c.) , (13) 

i=l,2 



Pll 



-7iPii - 2 (^1*^ "^'^Pis - c.c.[ 



hi = -72P22 - 2 (^2e"'''V23 - c.c) , 

PSl = -(1^31 -f 73i)p31 - 2^l6~"^'*(P33 - Pll) 

+ ^f]2e-"'^V2i, 

P32 = -(1^32 + 732)P32 " ^^2e~"'^*(p33 " P22) 

+ ^f]ie-"'^Vi2, 

P2I = -(it^21 +72l)P21 - 2^ie~"''V23 



(14) 
(15) 

I 

(16) 
) 
(17) 



"26 P31, 



(18) 



where c.c. implies the complex conjugate of the preceding 
term. It is useful to note here that these equations could 
be written more generally for a dense media in terms of 
the total local field. Within the linear response theory, 
and assuming the material under consideration is linear, 
we will take into account Lorentz-Lorenz correction after 
determining the dilute material response as usual |l5|. 

The relaxation rates of the off-diagonal elements of the 
density matrix are introduced as 27.^ = 7i + 7j . Formal 
solution for pi\(fy can be written as 



— lo*oi''2t 



i(t-i')e~''*"''+'''^+''''l*'di'. 



P2i(t) = 2^2e'-'^ / 93 
which leads to 



P31 = -(i^ + 73i)P3i - 2^i(/'33-Pii) 



^(i_i')e-i[(A-5)+72i]t'j^/^ 



p-i 



Here, we introduced a slow variable p3\ — /93ie"^^*, de- 
tuning of the probe beam A = (^31 — 1^1, and detuning of 
the driving beam 5 = 0^32 — v^. 

The effect of weak probe field on the system can be 
treated perturbatively. Carrying out standard linear re- 
sponse method, to the first order in the probe field am- 
plitude, we replace the inversion (/333 — pn) by its initial 
value which is taken to be —1, assuming only a small 
fraction of atoms are pumped out of their initial states. 
The integral can be evaluated by assuming ^31 doesn't 
change appreciably in time scale of I/721. We find 



P31 



-r^i 



[i(A - 5) + 721] 



(zA + 73i)[i(A - <5) + 721]+ I 1^2 P /4' 



(19) 



Positive frequency component of the complex induced 
electric dipole moment of the atom is given by -pi — 
d\^P3i, which is related to the complex atomic polariz- 
ability tensor a as Pi = aijEij. We adopt the summation 
convention, in which summation over a repeated index is 
implied. 



For the macroscopic polarization we have to take into 
account local field effects which lead to Clausius-Mossotti 
relation between polarizability and the susceptibility 
Xe- For small enough concentration TV of atoms Xe — 
Naeo holds. Using Pi = eox^J Eij = d\^p3i, we identify 
the complex electric susceptibility tensor for a gas of such 
three level atoms with concentration A^ to be [3, 13 



D 



2j3iheo D 

731 



nl 



Xe 



^ Nail 



N 

( 

3eo 



473i[i(A-(5)+72i] 
1 



(20) 



(21) 



Complex dielectric permittivity tensor can be similarly 
constructed via eij = €o{6ij + Xe)- We observe that 
this contributes to the complex permeability tensor of 
the system. It should be noted that for (5 = and for 
small N, we recover the well-known results for an elec- 
tromagnetically induced transparent system. Now, using 
the equation 

P21 = -(i^2i + 72i)p2i + ^r!*p3ie'(''^-'^i)*, (22) 
we deduce the relation 

i n% 



P21 



2 i(A -6)+ 721 



P31, 



(23) 



for the new variable p2i = P21 exp (i(i/i — 1^2)^- 

We can now calculate the induced magnetic dipole mo- 
ment of the atom using 



(/i) =Tr{pjl), 



(24) 



where fl = fisL/h is considered for the magnetic dipole 
operator by assuming the contribution from nuclear spin 
are negligible. The electronic spin part is for simplicity. 
As the lower levels are of opposite parity with the up- 
per level, the only non-vanishing contribution may arise 
if the lower levels are of the same parity. In this case we 
get (fl) = P21M12 + c.c, which gives 



iP) 



51il72^i2 exp (i(iyi - i^2)t) 
'4(iA + 73i)[i(A-(5) + 72i]+| ^^2 



-t- c.c. (25) 



In order to describe the atomic response to the mag- 
netic field component of the probe field, we let the in- 
duced magnetic dipole of the atom oscillate in phase with 
the probe beam. This is achieved when i^i — 1/2 — ^i^i- 
Setting aside the static field solution we consider the case 
of 1^2 = 2z^i. The other possibility j/2 = would be the 
case of a static electric field as the coupling field. This 
should be separately discussed as it is necessary to ex- 
amine Stark shifts of the levels and modify the present 
theory accordingly. The driving field is taken to be res- 
onant with the a;32 when the probe is resonant with the 



1^31 so that 5 — 2A which puts a constraint on the three 
level system as 0^32 = 2a;3i. 

This constraint is, however, a major obstacle in realiz- 
ing the predicted effects here at a realistic experimental 
setting as it is not straightforward to find a system with 
two states, which have a matrix element of fl between 
them and at the same time have energy difference in the 
optical range. This is mainly due to the fact that p. is an 
angular operator and the two states involved should have 
the same radial wavefunctions to give a non zero matrix 
element. One can imagine, some external magnetic field 
adjusting the separations to give the necessary energy 
conditions. However, for an atomic system to get split- 
tings in the optical regime, the external field would be 
impractically large. One can try to investigate systems 
in which fl is not an angular operator, such as molecular 
gases, and try to find optically separated states which 
have a magnetic dipole matrix element between them. 

As far as atomic gases are concerned, the best option 
seems to be to take two states which have the same L 
value but which are split due to L-S coupling to be the 
states |1) and |2) and try to get a third level of opposite 
parity to fulfill the energy condition. Another direction 
to proceed would be to consider all our discussion for an 
atomic system under high electric field. In that case, it 
will not be too hard to get to fields which give shifts on 
the order of optical frequencies, however one must care- 
fully do the preceding analysis again taking into account 
the effect of static electric field on all three states. 

We assume this condition is fulfilled with our hypo- 
thetical model atom and proceed by writing the prod- 
uct rii/zi2 explicitly, so that we can examine the direc- 
tional character of the magnetic response of the atom to 
the probe field. Electric dipole of the probe transition 
and the magnetic dipole of the lower levels are combined 
through a tensor product relation such that 



^1 i.fll2)i 



El 






-ij. 



where we introduce 



z.^^- = (1|m.|2)(3|4i|1 



(26) 



(27) 



The tensor v demonstrates the combined effect of elec- 
tric and magnetic field components of the optical field on 
the directional character of the magnetic response of the 
medium. 

To calculate the induced magnetic dipole moment ma- 
trix elements, it is convenient to consider angular mo- 
mentum basis in which we can also calculate the elements 
of electric dipole moment using the Wigner-Eckart theo- 
rem. Let us identify the states as 



|1) = \nj,m), 
|2) = \nj,m~l), 
|3) = \n',l + l,m- 1) 



(28) 



Using Lx = (i+ + L_)/2 and Ly = (L+ — L_)/2i matrix 
elements of the angular momentum are readily obtained 



in this basis as 



(1|L,|2) = -^{l + m){l-m + l), 



{l\Ly\2) = -^(l + m)il-m + l), 
(2|L.|1) = 0. 



(29) 



The matrix elements of the electric dipole operator can 
be conveniently calculated by expressing it as a spherical 
tensor operator of rank 1 so that its components become 



ez = Tf 



(1). 



1 



' 



.(1) ^(1) 



„ = -^(rv-ri 



1 I^tCi) ,r(i) 



ey = —\^ir{+T{-'J. (30) 

We use the Wigner-Eckart theorem in the form 



{nl3m3\T^l>\n'hmi) 






{nh\\T('-^\\n'h) 

V2I7TT 



Here, the first factor is the Clebsh-Gordan coefhcient 
where we choose a notation resembling its symmetric 
form in terms of Wigner-3j coefficients. The second fac- 
tor is the reduced matrix element which is independent 
of the orientation of the magnetic dipole characterized by 
the angular momentum projection quantum number m. 
In our case it is given by 



{n , I + l||er||n, I) 



drer'K,,+,ir)Rniir), (31) 



which is always non- vanishing, with Rni (r) being the ra- 
dial wavefunction. By the ?7i-selection rule (mi -|- 7712 = 
77i3 is required for non- vanishing matrix elements), we see 

that matrix elements of the T| g vanish. The sole non- 
vanishing matrix element of T_i determines the matrix 
elements of x and y components of the position operator 
which are found to be 



"31 



^'-'''+'^^'-^ + '^'n',l + l\\erMnl), 



y 2{2l + 2){2l + 3)^ 
4i = *c?3i: dl, = 0. 
Combining Eq|55| and Eq|21| we finally get 



A^B 



{Ti',l + l\\er\\nl){l-m + l) 



(32) 



(33) 



-1 


-i 0' 


i 


-1 









'{l + m){l -m + 2) 
{l + l)(2l + 3)^ 



as the matrix which determines the orientation of the in- 
duced dipole moment. It should be noted that the matrix 
v gives zero response to positively polarized EM waves 
in accordance with the dipole selection rules. For our 
particular set of levels, we need negatively polarized EM 
waves as they provide the photons with correct helicity to 



satisfy the angular momentum conservation in the probe 
photon emission and absorption processes between the 
states 1 and 3. For negatively polarized waves, i^ just 
reduces to a scalar. 

It is worth noticing that the structure of the tensor 
1/ resembles that of gyrotropic substances with both e 
and ^ are tensors such as pure ferromagnetic metals and 
semiconductors. These were argued to be most likely 
candidates to demonstrate left-handedness in the original 
paper by VeselagoQ- 

Further calculations require setting the polarizations of 
the coupling and the probe beams. To cause transitions 
between |2) and |3); the coupling beam polarization e^ 
has to have a component along the quantization direction 
z. So let us take the coupling beam to propagate in the 
X — y plane, and be linearly polarized along z. To cause 
transitions between states |1) and |3), the probe beam 
must have a polarization vector lying in the x — y plane. 
Let us take it to be propagating along the z axis with 
polarization lying in the x-y plane. 

Then, our general expression for the induced magnetic 
moment leads to 



fl^-f{Lu)E{x-iy), 



(34) 



where w now denotes the frequency of the probe beam 
and 



MBp*/„/ 

2n' 



jl^uj) = ^ni{n',l + l\\ef\\nl){l-m + l) 



{I + m){l - m + 2) 1 
""V (? + l)(2Z + 3)3 Z 

Z = 4(iA-H73i)[i(A-<5)+72i]+|^2|' 



(35) 



With this definition of 7(w) we can extend our result 
to macroscopic electromagnetics of a gas with concentra- 
tion N. In the spirit of Clausius-Mossoti equation [g, we 
define the magnetization per unit volume as 



M = N^{uj)(E + ^){x - ty) 
= Nj{Lu)il + ^)E{x-iy). 



(36) 



Now we recall the Fourier transform of the curl equation 
for electric field in Maxwell's equations. For a negatively 
polarized wave 



B 



1 



k X E{x — iy) 



= -E{x-iy). 



Combining equations (|36|l . (|37|) . we have 

M = -iN-f{uj)c{l + ^)B. 



(37) 



(38) 



Finally by using the definitions B = fj.o{H + M) and 
B = fiH, wc get 

1 



Mr(^) = 



1 + ZM07(^M1 + ^) 



(39) 
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FIG. 2: Frequency dependence of the relative (effective) di- 
electric permittivity tr and the relative magnetic permeability 
Hr of the dense gas of three level atoms with A'' — lO^^m"^, 
A ~ 589nm, 7 -- 10.06MHz, 7,,^ = 0.67, 'ygr/2Tv ~ lO^Hz, and 
^2 = 0.567. 





FIG. 3: Same with Fig|21but for the case of dilute gas with 
A*' — 10^ m~ . Here, the electric susceptibility and the mag- 
netic susceptibility are plotted as the relative permittivity and 
the permeability do not change appreciably from unity. 



with /.*,. = /i//io is the relative (effective) permeability. 
We shall see that combined effect of electric and magnetic 
field components of the optical fields, as well as local 
field effects lead to novel light propagation regimes in 
particular on the EIT resonance frequency. 



III. RESULTS AND DISCUSSIONS 

We consider a gas of ^^Na atoms with N — lO'^^'m^'^ 
to examine the case of dense media where the Lorentz- 
Lorenz local field corrections play significant role and 
N — lO^^m^"^, for the case of a dilute gas where the 
local field effects are weak. Our results are presented in 
Fig|21for the dense media and in Fig|31for the dilute gas. 

We see that both the relative dielectric permittivity 
Er = e/eo and the relative magnetic permeability /ii- = 
/i//io can become negative over a band of frequency ^ 
O.OOI7. This allows the propagation of light through oth- 
erwise opaque medium at high densities where the elec- 
tromagnetically induced transparency would not work. 
At resonance we find /ir(0) = —0.69 — iO.ll and £^(0) = 
— 1.86 + i0.12. It is natural to have transmission losses 
in our model, similar to other left-handed structures, as 
they are unavoidable due to the Kramers-Kronig rela- 
tions ensuring the causality in the system. On the other 
hand, theoretically it is not a trivial task to estimate the 
amount of losses ^^ [ITI I18| and to rigorously prove the 
causality in left-handed materials[T9, 20] . We can give 
a simple and rough estimate by simply taking into ac- 
count the imaginary part of the refractive index which 
gives that after several microns the optical field will be 
damped by ~ 33% due to linear absorption. At such 
length scales, our atomic system with the given densities 



may be found in Bose-Einstein condensed state due to 
the interatomic interactions. Multiple scattering of pho- 
tons as well as higher order many body correlations may 
contribute in addition to the local field correction. Such 
effects are argued to be about the same order with the 
local field correction [H EHH IH IH IH ll^ . It is an 
intriguing possibility that the present result of induced 
left-handedness could improve and benefit from contri- 
butions arising from the quantum correlations in a dense 
Bose-Einstein condensate or in a dense degenerate Fermi 
gas. In this paper, we will be content with limiting our- 
selves to classical gaseous media and hope to discuss the 
case of quantum gases elsewhere in detail. 

In the dilute gas limit, we recover the usual behavior 
of the electric susceptibility under electromagnetically in- 
duced transparency conditions. The transparency region 
is in the valley between the twin peaks in the imaginary 
part of the electric susceptibility where the peaks cor- 
respond to the two dressed absorption lines, the Autler- 
Townes doublet |23|. The magnetic susceptibility exhibits 
steep variation over a narrow band of frequencies in the 
vicinity of the resonance while in magnitude the rela- 
tive permeability remains close to unity for all A. When 
N ~ lO^*', similar results to those shown in Fig O are 
found where now Hr varies between 1.002 and 0.9985 over 
Ae (-7,7). 

In our numerical calculations, we estimate the dipole 
matrix clement from the spontaneous emission rate 7 ~ 
10.06MIIz using the relation ^31 — y^S^heo^^/STr'^. Here 
A is the wavelength of resonant probe transition which 
is A '^ 589nm. Typical values for 7ge = O.57 and 
jgr/2TT ~ lO^Hz are used. Rabi frequency associated 
with the driving field is chosen to be ^2 = O.567. 



IV. CONCLUSION 

In summary, we suggested a method for optical mod- 
ification of magnetic permeability using a three level 
scheme and derived the necessary conditions for its ap- 
plicability. We found that it is in principle possible 
to clectromagnetically induce Icft-handedness to a spa- 
tially homogeneous media. The major challenge we face 
is to have two levels separated at optical frequencies 
while having a non-vanishing magnetic dipole matrix cl- 
ement. Such level splittings require large external mag- 
netic fields or should be engineered by other means such 
as external electric fields or spin-orbital couplings. One 
may also consider molecular gases, or try to utilize ex- 
citonic energy levels in solid state heterostructures to 
engineer three level system fulfilling the energy condi- 
tion. The predicted effect is fundamentally based upon 
the Lorentz-Lorenz local field contribution in an elec- 
tromagnetically induced transparency medium of three 
level atoms with a non- vanishing dipole moment between 
lower levels. In dense medium limit, in which the medium 
becomes opaque normally with a negative dielectric con- 



stant, the presence of magnetic dipole gives rise to a nega- 
tive magnetic permeability so that the probe beam would 
still propagate within the otherwise optically thick dense 
medium for several microns before it is finally absorbed. 
It should be emphasized that the presented method 
is applicable to spatially homogeneous media and does 
not need any spatial periodicity which is unavoidable in 
metamaterials. In the dilute medium limit, the value 
of permeability do not change from unity appreciably, 
however, in this case we observed that it demonstrates 
steep changes over a small band of frequency. Such large 
gradient of permeability may affect the character of light 
propagation such as its group velocity and may serve an 
additional method to slow down or speed up the light. 
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